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SPECIFICATION 



IMAGE RECONSTRUCTION METHOD AND TOMOGRAPH 



Technical Field 

[0001] The present invention relates to a method for 
reconstructing a tomogram from a projection data in a 
computerized tomograph (hereinafter referred to as CT 
apparatus) using a fan beam or cone beam, more particularly to 
a method to reconstruct a tomogram from projection data having 

more than n [rad] of projection data angle using for 
back-projection and a CT apparatus to achieve the objective of 

the method thereof. 



Background Art 

[0002] An X-ray CT apparatus acquires projection data by 
receiving an X-ray radiated from an X-ray source and transmitted 
through an object being examined, with an X-ray detector 
disposed on the opposite side of the X-ray source interposing 
the object therebetween. Upon acquiring the projection data, 
the X-ray source and the X-ray detector disposed to be opposite 
to each other interposing the object therebetween are revolved 
around revolving axis, and the projection data in different 

• - • 

angles of rotation (phase) are collected. By reconstructing 
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this projection data, the creation of the internal image of the 
object is attained in a non-destructive manner. 
[0003] Such X-ray CT apparatuses are one of two kinds . One uses 
a single-array detector of which the detector elements are 
arranged one-dimensionally (line) , and another kind uses a 
multi-array detector of which the detector elements are 
arranged two-dimensionally . 

[0004] An imaging method in the simplest X-ray CT apparatus is 
the normal scanning method for creating an image by revolving 
an X-ray source and a detector around a revolving axis in a range 
of 2 7i , and the scanning range of the projection data acquired 
by this normal scanning method is 2 7c [rad] . 

[0005] On the other hand, since the above-mentioned fan beam 
or cone beam diverges centering around a central beam directed 
toward the revolving axis from the X-ray source, with regard 
to one beam directed toward the detector from the X-ray source, 
the same projection data (line integral) is measured twice 
during one revolution of the X-ray source and detector. Since 
such redundancy of data should be minimized to reduce the X-ray 
exposure, the imaging method for setting the scanning range at 

less than 2 tc is also adopted. As shown in Fig. 14 (a) and (b) , 
the beams are equivalent when the position of the X-ray source 

and the position of one detector element are switched, and by 
setting the maximum fan angle of the fan beam as 2 7m as shown 
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in Fig. 14 (c) , the projection data of all beams necessary for 
the image reconstruction can be measured at the point of the 

X-ray source moving by ^ +2 y m. This range is the minimum 
scanning range . 

[0006] However if the image reconstruction as in the case of 

scanning range is 2 n is carried out corresponding to the 
projection data obtainable from it +2 y m scanning range, the 
image gets distorted and the image quality deteriorates. This 
is due to the phase range of data possible to perform back 
proj ection being different with respect to each pixel . In other 
words, for example, as shown in Fig. 15, in pixel pi the data 

by which the phase range is more than tz is used for the image 
reconstruction centering around pixel pi, but in pixel p2 only 

the data by which the phase range is less than 7i is used 
centering around pixel p2 . 

[0007] This means that the redundancy of the projection data 
is different depending on the pixel. This is illustrated in 
a sinogram in Fig. 16 to indicate projection data, by 
representing the fan angle (an angle formed by the central beam 

and the respective beams) y on the horizontal axis and the 
revolving phase angle /? on the vertical axis. In other words, 
Fig. 16 is a sinogram showing a minimum complete data set, and 
the upper and lower triangle portions denoted with diagonal 

■m ' • ■ 

lines are the data being redundant to each other. 



As a method to solve the problem relating to the redundant 
data as mentioned above, it is suggested in Patent Document 1 
to assign weight to, for example, a predetermined region of the 
projection data. 

Patent Document: JP-A-2001-299738 
[0008] With regard to weighting function w for the fan beam it 
generally is required to fulfill formula (1) , and to fulfill 
formula (2) regarding weighting function w for the parallel 
beam . 

[Formulas 1] 

oo 

£ {w( p + 1%n, y) + w(n + p+2y + 2jtn-y)} =1 ( 1 ) 

71=0 

OO 

£{^(0 + 2tt«) + w(7t + j3 + 2m)}= 1 (2 ) 

However, while these weighting functions can be applied 
to the image reconstruction from a minimum complete data set 
(7t+2ym of scanning range) or the full-scanning data set (2 
71 of scanning range) , they are not applicable to the image 
reconstruction from a data set in which the scanning range is 
between these ranges. To solve this problem, a weighting 
function to be applied to data sets of intermediate-range is 



4 



suggested in Patent Document 1. Here, as shown in Fig. 17, by 

setting and using virtual fan angle r which is not dependent 
on the actual and physical maximum fan angle, the reconstruction 

from the projection data with a scanning range of 7c+2ym - 2 
7i is achieved. 

[0009] However, these weighting functions suggested in the past 

are not applicable to the scanning range of over 2 tc . Since 
the weighting function to apply for the scanning range changes 
to a different one at the point of the scanning range being over 

2 7i , the image quality such as noise quantity or artifact 
intensity will also be different between the result of the 

reconstruction from the range narrower than 2 7i and the result 
of reconstruction from the range wider than 2 7z . 

Also, on the virtual sinogram shown in Fig. 17, in the case 
that the different weight is assigned to the two triangles in 
the revolving phase direction, the configuration of the 
weighting function turns out to be a trapezium, triangle or the 
deformed non-linear shapes thereof, and the configuration gets 
closer to a triangle from a trapezium as the scanning range 

draws closer to 2 tc . This means that the region of which the 
weighting factor is less than 1 increases as the scanning range 

gets closer to 2 tz , and the data contribution ratio decreases 
significantly compared to the case that the scanning range is 

■ ■ • 

2 7t and the weighting factor of all the range thereof is 1 , which 
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can lead to a notable increase of noise (i.e. a decrease of SNR) . 
[0010] Another common problem is that the imaging noise 
decreases as the imaging data amount to be used for the 
reconstruction processing increases. In other words, the 
imaging noise decreases as the projection data width 
(projection data angle using for back-projection) increases. 
However, acquiring a wide phase range means the redundant 
imaging of the same place as shown in Fig. 18, which accompanies 
the decrease of the measurement through-put (spiral pitch, beam 
pitch and table feeding speed) . In this way, the decrease of 
the image noise and the reduction of imaging time are in 
trade-off relationship to each other, and the relationship 
between them can be inappropriate depending on the imaging 
purposes. 

[0011] On the other hand, in order to reduce the contradiction 
of data by motion movement and the deterioration of image 
quality caused by it, the arithmetic addition of the same data 
is performed. More specifically, for example in normal 
scanning, if there is no movement of the object being examined 
during one revolution of the X-ray source and detector, the 
projection data of imaging start-time phase (j3=0) coincides 
with the projection data of imaging end-time phase ( j3 =2 it ) . 
However, since it is impossible to completely remove the motions 
such as heart beats or blood flow, as shown in Fig. 19 (a) ~ (c) , 
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a discontinuity (data inconsistency) is generated in both 
projection data 51 and 52 which leads to a notable deterioration 
of image quantity such as streak artifacts 53 and 54. This 
discontinuity can be reduced by acquiring the identical data 
of imaging start time and imaging end time, and performing 
weighted addition between them. However, the assignment of a 
small weight causes a decrease of data contribution rate that 
leads to a decrease of SNR. In this way, the amount of image 
noise and corrective effect of discontinuity are in a trade-off 
relationship, and this relationship can be inappropriate in 
some situations . 

Disclosure of the Invention 
Problems 

[0012] As described above, with a conventional tomograph 
presents a difficulty to achieve the reconstruction in the 
projection data angle using for back-projection and also the 
problem that the discontinuity of the image is generated when 
the projection data angle using for back-projection is more than 

2 7i . It is also difficult to perform the correction of 
redundancy on the imaging data without a complicated procedure. 
A further problem with the conventional tomograph is the 
deterioration of image quality by discontinuity due to the 
movement of the object being examined, and the difficulty in 
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adjusting the trade-off relationship between reducing the 
amount of image noise and shortening the imaging time. 
[0013] The present invention was implemented in order to solve 
various problems that the above-mentioned conventional 
tomogram reconstruction method have, and its objective is to 
provide a reconstruction method using the weighting function 
applicable to the reconstruction of the projection data in the 
entire scanning range which enables the prevention of the 
lowering of the data contribution ratio and the acquisition of 
images with good SNR, and to reduce artifacts caused by the 
bio-movement. 

Means to Solve the Problems 

[0014] In order to achieve the above-mentioned objectives, the 
present invention provides an image reconstruction method by- 
creating a weighting function based on an arbitrary correction 
angle width (movement correction range) and projection data 
angle using for back-projection (width in the view direction 
to use for the reconstruction) taking a degree of bio-movement 
or redundancy into consideration, by using this weighting 
function . 

[0015] The image reconstruction method of the present invention 
comprises as follows. 

A tomogram reconstruction method which revolves the 
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radiation source and the detector disposed to be opposite to 
each other interposing a scanning object therebetween centering 
on a predetermined revolving axis, detects the penetrated 
radiation from the radiation source which filtered out the 
object being radiated, and creates a tomogram of the region of 
interest of the object from the detected projection data, 
including : 

a step for obtaining a weighting factor in compliance with 
the correction angle width and the projection data angle using 
for back-projection of the projection data; 

a step for obtaining the proj ection data to which the weight 
is assigned by carrying out the weighting process based on the 
weighting factor corresponding to said projection data; and 

a step for reconstructing the tomogram using the weighted 
projection data. 

[0016] The reconstruction method of the present invention can 
further include: 

a step for setting at least one of the following, a 
correction angle width and/or a projection data angle using for 
back-pro j ection ; and 

a step for setting the value of either the correction angle 
width or the projection data angle using for back-projection 
which was not set in the previous setting step, based on the 
value of the other one being set in the previous step. 
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The correction angle width and the projection data angle 
using for back-pro j ection are set, for example, as 0^ £ ^ (2F-1) 
(only £ ^2F-2 ceil<log2F> ) when the correction angle width is set 
as £ 7i and the projection data angle using for back-projection 

as 2F. Hereinafter s is referred to as the correction angle 
width index and F as the index of projection data angle using 

for back-pro j ection . 

[0017] The correction angle width is for eliminating the 
discontinuity due to motion movement being generated between 
a set of projection data and the following projection data to 
be measured and also for correcting data redundancy. The 

correction angle width index s is for representing the minimum 
guaranteed width of the slope portion of the weighting function . 
The correction angle width is set corresponding to the range 
of the region for correcting the data discontinuity 
(hereinafter referred to as the data discontinuity region) in, 
for example, the end portion of the projection data. Or, the 
correction angle width may be changed according to the amount 
of noise or motion artifact in reconstructed images. The 
correction angle width can also be increased or decreased by 
making it directly proportional to the projection data angle 
using for back-projection. 

The projection data angle using for back-projection (data 

» • * 

width) is determined considering the time-width (time 
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resolution) contributing to the data redundancy, SN, and images. 
It can also be set at any angle that is greater than data width 
of a minimum complete data set [two-times the value of the 
maximum of 7c + fan angle]. 

[0018] The weighting function being created based on the 

correction angle width index f and index F of projection data 
angle using for back-projection is set so that the weight of 
the data discontinuity region (weighting factor) turns out 
smaller than the weight of other region that is equivalent to 
said data discontinuity. For example, the first sub weighting 
function and the second sub weighting function being created 
by shifting the second sub weighting function by a predetermined 
phase can be added and normalized. 
[0019] To be more precise, it is obtained by: 

(1) First obtaining one weighting function (a first sub 
weighting function) which has the same size 
redundancy correcting region as the set bio-movement 
correcting range, which satisfies the previously 
mentioned formula (1) or (2) ; 

(2) making the weighting function which has the same form 
as the above-mentioned weighting function as a second 
sub weighting function by shifting it for a 
predetermined phase being determined by the scanning 
range ; and 
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(3) performing addition on the first sub weighting 
function and the second sub weighting function, and 
normalizing it (performing the multiplication on the 
gain for sub-weight function) . 
As an example of the sub weighting function, a shape of 

trapezium with tz - e tz as an upper hem and tz + e tz as a bottom 
can be cited. 

[0020] The projection data applying the weighting function 
obtained in this way has no distortion of images caused by 
redundancy, due to the discontinuity among data being corrected 
according to the demand of users indicated with the correction 

angle width index. Also if the scanning range is more than tz , 
it can be applied to the projection data with any projection 
data angle using for back-projection whether it is less than 
2 tz or more than 2 tz , so the problem with the discontinuity of 
which existed at 2 tz of scanning range is solved. 
[0021] Also, the reconstruction method for the tomogram in the 
present invention further includes a step for implementing the 
process to rearrange the fan beam radiated from the radiation 

source to the parallel beam, and the weighting functions w( 6 ) 
when the revolving phase (view phase) at the time of detecting 

the pro j ection data is set as 9 and using N which can be obtained 
from the correction angle width e tz [rad] and 2 (N_1> ^ F- f /2<2 N 
(N is an integer of more than 0) are: 
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[Formula 2A] 



w( 6 )=0 



if [ 0 < Po 7c ] 



W( 0 ) = (p 7 7C+ 0 )W1/ (f 7C) 



if [Po 7t ^ 6 <p x 7z , e >0] 



w( 0 )=0 



if [PO 7C ^ 0 <P X 7C , £ =0] 



w( 0 ) =W1*V2*2/ £ 



if [PiTC ^ 0 <p 2 7C , £ >0 f V1 = 0] 



w 



( 0 ) = ( ( 0 -Pi 7i ) * (Wl*4/ £ ) /2 7C ) +W1*V2*2/ £ 



if [Pi7C ^ 0 <P 2 7i y e >0 ,V1=£0] 



w ( 0 ) =W1 



if [PiTC ^ 0 <P 2 tu , £ =0] 



w ( 0 ) = ( ( 0 -P 3 ^ ) *wi/ ( e 71 ) ) +W2 



if [P 2 tc ^ 0 <P 3 tc , £ >0] 



w ( 0 ) =W2 



if [P 3 7C ^ 0 <p 4 7C ] 



W ( 0 ) = ( (P 4 71 - 0 ) *W1/ ( € 7C ) ) +W2 



if [P 4 7C ^ 0 <P 5 7C # £ >Q] 



w ( 0 ) =W1*V2*2/ e 



if [P 5 7C ^ 0 <P 6 7C , £ >0,V1 = 0] 



w ( 0 ) = ( (P6 7i - 6 ) * (Wl*4/ f ) /2 7T ) +W1*V2*2/ £ 



if [P 5 tc ^ 0 <P 6 7c , £ >0, VI =£0] 
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w ( 0 ) =W1 



w( 0 ) = (P 7 7t-d )W1/ (en) 



w( 0 )=0 



w( 0 )=0 



if [P 5 7C ^ 0 <P 6 7C f £ =0] 



if [P 6 7t ^ 0 <P 7 7C , £ >0] 



if [P 6 n S 0 <P 7 n , £ =0] 



if [P 7 7C ^ 0 ] . 



Here, the respective parameters in the above formulas 
determined by the following respective formulas : 
[Formula 2B] 



Vl= £ -F + 2 



(N-1) 



if [ £ -F+2 (N_1) >0] 



V2- £ /2-V1 
M=2 N 
W1=1/2 N 
W2 = 1/2 (N ~ 1} 
W2=2*W1 

W2= (2* (M-F) + £ ) *W1/ £ +W1 
AA=-F 

BB=-F+ £ 
CC=M-F 

DD=M-F+ £ 



if [ £ ^0] 



if [ £ >0 , F<M] 
if [ £ >0, M^F] 



EE=F-M- £ 
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FF=F-M 

GG=F- £ 

HH=F 

Po=AA 



Pi=BB 


if 


[F<M/2+ £ /2] 


Pi=EE 


if 


[M/2+ £ /2^F<M/2+ £ ] 


Pi=BB 


if 


[M/2+ £ ^F] 


P 2 =BB 


if 


[M/2 + £ /2^F<M/2+ £ ] 


P 2 =EE 


if 


[M/2 + £ ^F<M+ £ /2] 


P 2 =CC 


if 


[M+ £ /2^F] 


P 3 =FF 


if 


[M/2+ £ /2^F<M] 


p 3 =cc 


if 


[M = F<M+ £ /2] 


P 3 =EE 


if 


[M+ £ /2^F] 


P 4 =CC 


if 


[M/2 + £ /2^F<M] 


P 4 =FF 


if 


[M^F<M+ £ /2] 


P 4 =DD 


if 


[M+ £ /2 = F ] 


P 5 =GG 


if 


[M/2+ £ /2^F<M/2+ £ ] 


P 5 — DD 


if 


[M/2 + £ ^F<M+ £ /2] 


P 5 =FF 


if 


[M+ £ /2^F] 


P 6 =GG 


if 


[F<M/2+ £ /2] 


P 6 =DD 


if 


[M/2 + £ /2^F<M/2+ £ ] 


P 6 =GG 


if 


[M/2 + £ ^F] 


P 7 =HH . 
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[0022] The above-mentioned weighting functions can be expressed 
with formula (3) using the following formulas (4) - (8) : 



[Formulas 3 ] 



r 



Wp(6) = G 



Ws 



V 



e_ 



+ Ws 



J 



e 

2% 



J- 



(3) 



0 



^ (T? + e/2) 
2 



e/2 



— \H 

V 2 J 



J M 2 
otherwise, 



(4) 



0c - 5 



(5) 



_27j + e-F 

y c2 « 



(6) 



77 = 2"" 1 



(7) 



G = 2' N . (8) 

Here, W s is the sub weight for the parallel beam, 6 is the 
view phase, 9 cl and 0 C 2 are the central view phase of the sub 
weight, v is the sub weight reference width, G is the gain for 
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sub-weight function, N is an integer of more than 0 being 2 



(N-l) 



= F— £ /2<2 N . Correction angle width index e and/or back 
projection width index F are inputted by the user, and the other 
is determined automatically. 

[0023] Or, it also is possible to make formula (3) a non-linear 
weighting function by correcting the sub weight of it to be 
non-linear. The non-linear weighting function is expressed 
with formula (3') using the above-mentioned formulas (4) ~ (9) . 



[Formulas 4] 



Wp(6) = G 





( ( 


< NL 


Ws 







e 



271 



-G cl ,-n,e 







( ( 


c 


+ NL 


Ws 






K v 



6 



2% 



J 



J 



(3') 



NL(w) = 3w 2 - 2w 



(9) 



[0024] As for a fan beam, the weighting functions w (6 , y ) for 
the reconstruction of the fan beam in the case of setting the 

projection phase of the fan beam as 9 and the fan angle as y , 
using N which can be obtained from the correction angle width 

£ 7i [rad] and 2 <N_1) F- s /2<2 N (N is an integer of more than 0) 
are : 



[Formula 5A] 
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w( 0 )=0 



if [ 0 < Po 7i ] 



W( 0 ) = (P 7 K+9 )Wl/( f 7C) 



if [POTC ^ 0 <p x 7C f £ >0] 



W( 0 )=0 



if [POTC ^ 0 <P 1 7C , £=0] 



w 



( 0 ) =W1*V2*2/ £ 



if [PiTC ^ 0 <P 2 tc f £ >0,V1 = 0] 



w 



( 0 ) = ( ( 0 -P x 7t ) * (Wl*4/ £ ) /2 ) +Wl*V2*2/ £ 



if [PiTC ^ 0 <P 2 7C, £>0,V1=£0] 



w ( 0 ) =W1 



if [PiTC ^ 0 <P 2 7C f £ =0] 



w 



( 0 ) = ( ( 0 -P 3 7t ) *W1 / ( £ 71) ) +W2 



if [P 2 71 ^ 0 <P 3 7C , £ >0] 



w ( 0 ) =W2 



if [P 3 71 ^ 0 <P 4 7C ] 



w ( 0 ) = ( (P 4 7t - 0 ) *wi/ ( e 7i)) +W2 



if [P 4 7r ^ 0 <P 5 tc , £ >0] 



w( 0 )=W1*V2*2/ £ 



if [P 5 7C ^ 0 <P 6 n , £ >0 , V1=0] 



w 



( 0 ) = ( (P6 7i - 6 ) * (Wl*4/ £ ) /2n) +W1*V2*2/ £ 



if [P 5 7t ^ 0 <P 6 tt , £ >0,V1=£0] 



w ( 0 ) =W1 



if [P 5 7C ^ 0 <P 6 7C , £ =0] 
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w 



( 0 ) = (P? rc- 9 ) Wl/ ( s ti: ) 



if [P 6 7C ^ 0 <P 7 71 , £ >0] 



w( 0 )=0 



if [P 6 7C ^ 0 <P 7 71 , £ =0] 



w( 6 )=0 



if [P 7 tc ^ 0 ] . 



Here, the respective parameters in the above formulas are 
determined by the following respective formulas : 
[Formula 5B] 

« 

Vl= £ -F+2 (N ~ 1} if [ £ -F+2 (N_1) >0] 

V2- £ /2-V1 
M=2 N 

W1=1/2 N 

W2 = 1/2 <N ~ 1) if [ £ ^0] 

W2=2*W1 if [ £ >0, F<M] 

W2= (2* (M-F) + £ ) *W1/ £ +W1 if [ £ >0 , M^F] 



AA=-F 



BB=-F+ £ 



CC=M-F 



DD=M-F+ £ 



EE=F-M- £ 



FF=F-M 
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GG=F- £ 






HH=F 






Po=AA 






Pi=BB 


if 


[F<M/2+ £ /2] 


Pi— EE 


if 


[M/2+ £ /2^F<M/2+ £ ] 


Pi=BB 


if 


[M/2+ £ ^F] 


P 2 =BB 


if 


[M/2+ £ /2^F<M/2+ £ ] 


P 2 =EE 


if 


[M/2 + £ ^F<M+ £ /2] 


p 2 =cc 


if 


[M+ £ /2^F] 


P 3 -FF 


if 


[M/2+ £ /2^F<M] 


p 3 =cc 


if 


[M^F<M+ £ /2] 


P 3 =EE 


if 


[M+ £ /2 ^F] 


P 4 =CC 


if 


[M/2+ £ /2^F<M] 


P 4 =FF 


if 


[M^F<M+ £ /2] 


P 4 =DD 


if 


[M+ £ /2^F] 


P 5 =GG 


if 


[M/2 + £ /2^F<M/2+ £ ] 


P 5 =DD 


if 


[M/2+ £ ^F<M+ £ /2] 


P 5 =FF 


if 


[M+ £ /2^F] 


P 6 =GG 


if 


[F<M/2+ £ /2] 


P 6 =DD 


if 


[M/2 + £ /2^F<M/2 + £ ] 


P 6 =GG 


if 


[M/2+ £ ^F] 


P 7 =HH. 







[0025] Weighting function Wf ( 6 , y ) in relation to the fan beam 
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can be expressed with formula (10) by using the same presenting 

method as weighting function Wp ( 9 ) regarding the 
above-mentioned parallel beam, and also with formula (10') 
using the non-linear sub weight. Either one of these may be 
used with regard to the fan beam. 



[Formulas 6] 



wf(e,Y) = GAws 



(0-y 



v 



2% 



\ 



r 



+ Ws 



J 



V 



0-y 

2% 



\ 



(10) 



Wf(0 9 y) = G-\NL 





( ( 


[nl 


Ws 




I V 



2tc 



-O cl ,rj,e 







( ( 




+ NL 


Ws 


)) 







fe-y ^ 



2% 



v 



(10') 



[0026] Also, in the reconstruction method for tomograms in the 
present invention, the projection data can be data detected as 
an object being shifted in the revolving axis direction along 
with the revolution of the radiation source and the detector. 
In such a case, a step of interpolating the projection data and 
creating the projection data of the orthogonal side to the 
revolving axis (called "2-dimensional interval projection 
processing") is included. 

[0027] A tomograph of the present invention includes: 



a radiation source and a detector disposed to be opposite 
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to each other interposing an object therebetween; 

a reconstruction means for creating a tomographic image 
of the region of interest of the object from the projection data 
detected by the detector; and 

an imaging control means for controlling the radiation 
source, detector and reconstruction means, 

wherein the reconstruction means is mounted with the 
above-mentioned reconstruction method of tomograms. 

In the tomograph of the present invention, the imaging 
control means can change the projection data angle using for 
back-projection according to the imaging region. For example, 
it can improve SNR for imaging by widening the projection data 
angle using for back-projection, or can improve the time 
resolution for imaging by narrowing the projection data angle 
using for back-projection. Also, the tomograph of the present 
invention is provided with a meaning for moving an object in 
the revolving axis direction, wherein the imaging control means 
can change the correction angle width and/or the projection data 
angle using for back-projection to use for the reconstruction 
according to the movement velocity of the object in the 
revolving axis direction. 

[0028] In the tomograph of the present invention, the detector 
may be either a single-array detector or a multi-array detector . 
When it is the multi-array detector, the reconstruction means 
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may use the same weighting factor with regard to the respective 
array of the detector, or may use a different weighting factor 
from other arrays with regard to at least one of the rows of 
the detector. 

Moreover, the tomograph of the present invention 
preferably comprises an input means for accepting the 
information from the users relating to the correction angle 
width and the projection data angle using for back-projection. 

Advantageous Effect of the Invention 

[0029] According to the present invention, it is possible to 

i 

perform a reconstruction using the weighting function which can 
correspond to the back projection with either fan beam or 
parallel beam, and is also applicable to the consecutive phase 
width from narrow projection data angle using for 
back-projection having raised contribution rate to wide 
projection data angle using for back-projection having data 
redundancy, and to attain high-definition tomograms without the 
problems of data discontinuity or noise. 

[0030] Concretely, the weighting function in the present 
invention does not form a rectangle due to the minimum 
correction angle width being maintained by the correction angle 

width index e being set at more than 0 ( e >0) . Taking a 
weighting function shape for the parallel beam for example, the 
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conventional idea would be that the weighting function shape 
turns out to be a rectangle as seen in Fig. 11 (b) when index 
F of projection data angle using for back-projection =1.0, but 
with the same index F of projection data angle using for 
back-pro j ection = 1.0 by setting the correction angle width 

index £ at 0 . 2 the weighting function shape of the present 
invention turns out to be a two-tier trapezium as shown in Fig. 
11 (a) , and the weight slope is being maintained. This 
indicates that the influence caused by the movement can be 
reduced by using the weighting function of the present invention 
since the data discontinuity upon the phase being 0 and the phase 

being 2 tz is restrained by the data of phased which is the 
opposite data of the data discontinuity. Making the weighting 
function nonlinear enables the attainment of better results. 

Best Mode for Carrying Out the Invention 

[0031] The embodiment of the present invention will be described 
below referring to the attached drawings. 

Fig. 1 is a diagram showing an overview of a CT apparatus 
for implementing the reconstruction method of to tomograms of 
the present invention. This CT apparatus mainly comprises 
scanner 40 , operation unit 50 , and table 60 for laying and moving 
an object. Scanner 40 comprises units such as central 
controller 400, X-ray controller 401, high-voltage generator 
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402 , high-voltage switching unit 403, X-ray generator 404, 
X-ray detector 405, preamplifier 406, scanner controller 407, 
driver 408, collimator controller 409, table controller 410, 
and table-feeding measurement unit 411. 

[0032] Operation unit 50 comprises input/output unit 51 
including units such as the display unit, input unit and storage 
unit and calculation unit 52 including units such as 
reconstruction calculator and image processor. The input unit 
comprises a device such as mouse or keyboard, and is for 
inputting the measurement/reconstruction parameter such as the 
table feeding speed information or reconstruction position. 
In the present invention it is possible for a user to set 

correction angle width index s and index F of projection data 
angle using for back-projection as a necessary parameter to 
obtain the weighting function to use . for the image 
reconstruction of the projection data which is obtained 
additionally, and these index are also inputted from input unit. 

A detailed description of correction angle width index f and 
index F of projection data angle using for back-projection will 
be presented later . The storage unit is for storing information 
being inputted from input unit or the processing results in 
calculation unit 52. The display unit is for displaying the 
various types of data such as above-mentioned information or 
reconstructed images . The reconstruction calculator is for 
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processing the data obtained from the X-ray detector, and the 
image processor is for implementing various processing on 
images such as a reconstructed image and displaying them to the 
display unit. 

[0033] Central processing controller 400 is for transmitting 
the necessary controlling signals for imaging to X-ray 
controller 401, table controller 410 and scanner controller 407 
based on the imaging condition (such as table feeding speed, 
tube-current, tube-voltage and slice position) being inputted 
from the input unit of operation unit 50 or the reconstruction 
parameter (such as the region of interest, size of reconstructed 
images, projection data angle using for back-projection and 
reconstruction filter function) , and starting the imaging upon 
receiving the starting-signal. Upon starting the imaging, the 
controlling signal is transmitted to high-voltage generator 402 
by X-ray controller 401, high-voltage is applied to X-ray 
generator 404 via high-voltage switching unit 403, the X-ray 
outputted from X-ray generator 404 is irradiated to an object 
being examined, and the transmitted X-ray is inputted to X-ray 
detector 405. Simultaneously, a controlling signal is 
transmitted from scanner controller 407 to driver 408 , and X-ray 
generator 404, X-ray detector 405 and preamplifier 406 are 
controlled to revolve around the object. When executing 
circular orbit scanner (normal scanning) 41 as shown in Fig. 
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2(a), table 60 where the obj ect is laid down , comes to rest during 
the revolution at the command of table controller 410. Also 
when executing helical orbit scanning 42 as shown in Fig. 2 (b) , 
the location of table60 is shifted in the revolving axis 
direction of units such as X-ray generator 404 with a 
predetermined pitch. In the helical scanning method, when the 
distance of the table feeding relative to the scanner during 

one revolution of the scanner is set as distance Ax, the ratio 
of the detecting element in the revolting axis direction width 

D (D/Ax) is defined as the helical pitch, and the ratio of the 
detector to the entire length of the revolving axis direction 
is defined as the beam pitch. The higher these pitches are, 
indicates that the imaging time for the same range can be 
shortened and imaging performance can be improved. Generally 
the helical pitch is used to a degree of "2" which can cover 
almost the entire imaging region in consideration of the 
opposite data. 

[0034] The X-ray outputted from X-ray generator 404 gets 
restricted in its irradiation area by collimator 412 under the 
control of collimator controller 409, transmits through the 
object, being absorbed (attenuated) in the respective tissues 
in the body of the object, and gets detected by X-ray detector 
405. The X-ray being detected by X-ray detector 405 is 
converted into a current at that point, amplified by 
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f 



preamplifier 406, and inputted to calculator 52 of operation 
unit 50 as a proj ection data signal . The proj ection data signal 
inputted to calculator 52 receives the image reconstruction 
processing by the reconstruction calculator in calculator 52. 
This reconstructed image is stored in the storage unit of 
input/output unit 51, and displayed on the display unit as a 
CT image . 

[0035] X-ray detector 405 may be either single-array detector 
11 as shown in Fig. 3 (a) in which the detecting elements are 
arranged in one dimension, or multi-array detector 12 as shown 
in Fig. 3 (b) in which the single-array detector is arranged 
in a plurality of rows in the revolving axis direction (the 
direction of the arrow in the diagram) . Additionally, while 
detector 11 is shown in a straight line in Fig. 3 (a) , generally 
the one with the detecting elements arranged in a circular arc 
shape is used so that the distance from X-ray source 10 to the 
respective detector elements or the angle between the adjacent 
X-ray beams will be equal. In the single-array detector, the 
X-ray beam is orthogonal to the revolving axis. In the 
multi-array detector, a wide range can be imaged in one imaging 
compare to the single-array detector, but as it moves away from 
the mid-lane (central row) in the- revolving axis direction the 
X-ray beam forms a tilt angle (cone angle) . 

[0036] The processes to be carried out in the reconstruction 
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calculator are, for example, the rearrangement process for 
coordinating the fan beam projection data obtained by the fan 
beam with the parallel beam projection data in a parallel form, 
the data correction process for calculating the weighting 
function to apply to the fan beam projection data or parallel 
beam projection data and for applying the weighting function 
to these projection data, the filter correction process for 
generating the filter processing parallel beam projection data 
by superimposing the reconstruction filter to the parallel 
projection data, the back projection process for implementing 
the back projection of the filter processing parallel beam 
projection data to the back projection region relating to the 
region of interest, and the data interpolation process for 
creating the data like circular orbit by carrying out the data 
interpolation to the data obtained upon the helical orbit 
scanning. Also in the case that detector 405 is the multi-array 
detector as shown in Fig. 3 (b) , the cone angle correction 
process for performing the multiplication on the coefficient 
depending on the tilt angle of the radiation relating to the 
respective row of projection data is carried out. 
[0037] Next, the tomogram reconstruction method in the X-ray 
CT apparatus of the above-mentioned configuration will be 
described. The procedure of the reconstruction method is 
illustrated in Fig. 4 (a) . The imaging method may be either 
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one of the circular orbit scanning or the helical scanning, and 
the case with the circular orbit scanning will be described 
here . 

[003 8] As mentioned above, the projection data can be obtained 
by X-ray generator 404, X-ray detector 405 and priamplifier 406 
revolving around an object while an X-ray outputted from X-ray 
generator 404 being irradiated to the object and the transmitted 
X-ray is imputed to X-ray detector 405 (step 101) . The scanning 
range is 2F7C being determined by index F of projection data 
angle using for back-projection inputted from the input unit, 
and in the case for the parallel beam it is more than n . in 
other words the numeric value that satisfies F^0.5 as the index 
F of pro j ection data angle using for back-pro j ection is inputted . 
As for the circular orbit scanning, the user determines the 
appropriate value in consideration of the exposure quantity of 
the X-ray, time resolution or noise. 

[0039] First, the obtained fan beam projection data is 
rearranged to the parallel beam projection data (step 102) . 
This rearrangement process is the process of combining the data 
from different phases and fan angles from the fan beam being 
irradiated radially from the revolving axis direction as shown 
in Fig. 5 (a) , and converting it into a parallel beam that is 

parallel from the revolving axis direction as shown in Fig. 5 

■ 

(b) . SI and S2 in Fig. 5 indicate the position of the radiation 
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source and the detector. While the address on the detector of 
the X-ray beam passing through the reconstruction pixels in the 
respective projection phase is calculated in the back 
projection process that is heavily used as the reconstruction 
method, by using the parallel beam, the inverse trigonometric 
function which is time consuming and high in calculation load 
and the distance calculation between the X-ray source and the 
reconstruction pixel become unnecessary, and an inverse 
trigonometric function can be replaced with a product-sum 
operation. Therefore in spite of the increase of calculation 
time for the rearrangement process, it has the advantage of 
speeding up the time required for the reconstruction process . 
[0040] Next as the data correction process, a correction is 
carried out for eliminating image distortion due to the 
bio-movement correction of the projection data or redundancy 
(step 103) . In the data correction process, weighting function 
w(0) for applying to the parallel beam is obtained based on 
correction angle width index s and index F of projection data 
angle using for back-projection inputted from the input unit. 

Correction angle width index e is set by the user for the 
purpose of eliminating data discontinuity by bio-movement in 
consideration of the degree of bio-movement in compliance with 
conditions such as the imaging region of the object, and is set 
to be 0^ f ^ (2F-1) which is the range not surpassing [scanning 
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width (2F7c)-7c] when the correction angle width is set as e 

71 . 

[0041] The weighting function w ( 9 ) (9 is the pro j ection phase) 
satisfies the following formula (1) (in the case of a fan beam) 
or formula (2) (in the case of a parallel beam) with regard to 
the correction region of the projection data determined by 
correction angle width £ n as shown in Fig. 6, 

[Formulas 7 ] 



and is obtained by acquiring the sub weighting function W s ( 9 ) 
of the trapezium shape with n + a n as the bottom and k - s n as 
the upper hem, adding the function of which this sub weight 

function W s ( 9 ) being shifted to (2F 7i - ( n + e n ) ) by a 
predetermined phase being determined by scanning width 2F7C in 
relation to the center of data width, and implementing the 
normalization. Normalization is the process of making the 
average weight of the respective phase equal 1. 
[0042] In concrete terms , it is calculated by the formulas below 



£{m>( P + 27r«,7) + w{n + p + 2y + Inn ,-/)}= 1 



(1) 




(2) 
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using N which satisfies 2 <N 1} ^ F- £ /2<2 N (N is an integer of 
more than 0) . 



[Formula 8A] 



w( 9 )=0 



if [9 < Potc ] 



w 



( 0 ) = (p 7 it + 9 ) wi/ ( e 7t ) 



if [Pore ^ 9 <p x 7i , e >0] 



w( 9 )=0 



if [PO 7C ^ 0 <P X TT , € =0] 



w ( 9 ) =W1*V2*2/ £ 



if [PiTc ^ 9 <P 2 7t , e >0,V1=0] 



w ( 9 ) = ( ( 9 -Pi 7i ) * (Wl*4/ £ ) /2 ) +W1*V2*2/ £ 



if [Pi 7C ^ 8 <p 2 7C , £>0,V1=£0] 



w ( 9 ) =wi 



if [P 1 7l ^ 9 <P 2 7Z , a =0] 



w ( 9 ) = ( ( 9 -P 3 n ) *W1/ ( £ 7i ) ) +W2 



if [P 2 7C ^ 0 <p 3 7C , £ >0] 



w ( 9 ) =W2 



if [P 3 ^ = 9 <P 4 7C ] 



w ( 0 ) = ( (P 4 7Z - 6 ) *W1/ ( £ 71 ) ) +W2 



if [P 4 7t ^ 0 <P 5 7C , £ >0] 



w ( 9 ) =W1*V2*2/ £ 



if [P 5 7l^t 9 <P 6 7C , £ >0,V1 = 0] 
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w( 0 ) = ( (P6 7i- 0 ) * (Wl*4/ s ) /2 n ) +W1*V2*2/ £ 



w ( 0 ) =W1 



w ( 0 ) = (P 7 tc - 0 )W1/ ( e 7i) 



w( 0 ) =0 



w( 0 )=0 



if [P 5 7C ^ 0 <P 6 7C f € >0, VI ^=0 ] 



if [P 5 ?r ^ 0 <P 6 tc , £ =0] 



if [P 6 71 ^ 0 <P 7 71 , £ >Q] 



if [P 6 7T ^ 0 <p 7 7c , e =0] 



if [P 7 7T ^ 0 ] 



Here, the respective parameters in the above formulas are 
determined by the following respective formulas: 



[Formula 8B] 



Vl= a -F+2 



(N-l) 



if [ £ -F + 2 (N_1> >0] 



V2- £ /2-V1 
M=2 N 

W1=1/2 N 
W2 = 1/2 (N_1) 
W2=2*W1 

W2= (2* (M-F) + £ ) *W1/ £ +W1 



if [ £ ^0] 



if [ f >0 , F<M] 
if [ f >0 , M^F] 



AA=-F 



BB=-F+ £ 
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COM-F 
DD=M-F+ £ 

EE=F-M- £ 
FF=F-M 

GG=F- £ 

HH==F 

Po=AA 



Pi-BB 


if 


[F<M/2+ £ /2] 


Pi=EE 


if 


[M/2+ £ /2^F<M/2 + £ ] 


Pi=BB 


if 


[M/2+ £ ^F] 


P 2 =BB 


if 


[M/2+ £ /2^F<M/2 + £ ] 


P 2 =EE 


if 


[M/2 + £ ^F<M+ £ /2 ] 


P 2 =CC 


if 


[M+ £ /2^F] 


P 3 =FF 


if 


[M/2+ £ /2^F<M] 


p 3 =cc 


if 


[M^F<M+ £ /2] 


P 3 =EE 


if 


[M+ £ /2^F] 


p 4 =cc 


if 


[M/2 + £ /2^F<M] 


P 4 =FF 


if 


[M^F<M+ £ /2] 


P 4 =DD 


if 


[M+ £ /2^F] 


P 5 =GG 


if 


[M/2+ £ /2^F<M/2+ £ ] 


P 5 =DD 


if 


[M/2 + £ ^F<M+ £ /2] 


P 5 =FF 


if 


[M+ £ /2^F] 


P 6 =GG 

* 


if 


[F<M/2+ £ /2] 


P 6 =DD 


if 


[M/2+ £ /2^F<M/2+ £ ] 
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P 6 =GG 
P 7 =HH . 



if [M/2+ £ ^F] 



[0043] This weight function W P for the parallel beam can be 
expressed by formula (3) using the following formulas (4) ~ (8). 



[Formulas 9] 



Wp(6) = G\Ws 



re 



v 



2n 



-9 A ,%e. \ + Ws 



e 



v 



2% 



J 



(3) 



0 



e/2 



■q + e/2 



e 



if |ga (T?+e/2) 

otherwise, 



(4) 



0 



cl 



__2r\ + e — F 



(5) 



_ 2t] + e-F 

y c2 r 



(6) 



77 = 2 



(7) 



G = 2 



(8) 
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Here, W s is the sub weight for the parallel beam, 0 is the 
view phase, 6 cl and 8 C 2 are the central view phase of the sub 
weight, v is the sub weight reference width, G is the gain for 
sub-weight function, N is an integer of more than 0 to be 2 (N_1> 

= F- £ /2<2 N . Correction angle width index f and/or back 
projection width index F is inputted by the user, and the other 
is determined automatically. 

[0044] Or, it is also possible to correct formula (3) in such 
a way that the sub weight turns out to be non-linear, and make 
it a non-linear weighting function. In other words the 
non-linear weighting function can be expressed as formula (3') 
using the above-described formulas (4) ~ (8) and formula (9) 
described below. 



[Formulas 10] 





* 


( ( 


Wp(0) = G< 


NL 


Ws 









e 



2k 







( ( 




+ NL 


Ws 


J) 




\ v 



2k 



(3') 



NL(w) = 3w 2 -2w 3 (9) 



[0045] In concrete terms, the weighting function in the case 
when the correction angle width index is s =0 . 2 and when e =0 . 8 
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is illustrated in Figs. 7 and 8. In Figs. 7 and 8, each of (a) , 
(b) and (c) shows a case when index F projection data angle using 
for back-projection is 0.9, 1.0 and 1.1 respectively. As 
indicated in the drawings, this weighting function can be 

applied to the projection data of more than n and of arbitrary- 
projection data angle using for back-projection by setting 
index F of projection data angle using for back-projection. 

Also by setting correction angle width index f , it is possible 
to obtain definite corrective effect which enables correction 
of the discontinuity in the end portion of the data, even in 
the case of the projection data angle using for back-projection 

having any value such as 1 . 3 7t or 2 tz . Furthermore, this 
weighting function can use the same weighting function 
corresponding to all the channel positions of the parallel beam 
(the position corresponding to the fan angle position of the 
fan beam) , and store the weight with a small amount of memory. 
Also , the method according to the present embodiment can include 
redundant data in CT reconstruction without using the weight 
of which the phase width is limited. In the case of correction 

angle width index e = 0 , it is possible to obtain the projection 
image placing as much value on the reduction of noise as in the 
past . 

[0046] Next, with regard to the parallel beam projection data 
applying the weighting function being obtained in this way, the 
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reconstruction filtering process is carried out superposing the 
reconstruction filtering function in the channel direction 
(step 104), after that the back projection process is 
implemented and the CT image is obtained (step 105) . As for 
the reconstruction filtering and the back projection process, 
heretofore known methods can be used. 

[0047] While CT apparatus using a single-array detector was 
described in the above-mentioned embodiment, in the CT 
apparatus using the multi-array detector, since the projection 
data relating to the detection elements of the plurality of 
arrays placed in the revolving axis direction can be obtained, 
each of the above-mentioned rearrangement process (step 102) 
and correction process using the weighting function (step 103) 
are implemented corresponding to the projection data with 
respect to each row. Then by performing the multiplication of 
the cosine of the several cone angles corresponding to each row 
on the parallel beam projection data, carrying out the 
reconstruction filtering process which superimposes the 
reconstruction filtering function in the channel direction of 
the each row of parallel beam, and further implementing the 
3-dimensional back projection process, CT images similar to CT 
apparatus using the single-array detector can be obtained. 
Here the weighting function applied to the projection data of 
each row may be the same or different. The correction angle 
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width index can be set with respect to the movement of the imaging 
region. For example, in the case of imaging a plurality of 

regions by normal scanning, projection data angle 2F7C using 
for back-projection is widened in the region giving greater 
emphasis on noise amount and narrowed in the region giving 
greater emphasis on the time resolution. Also since data 
redundancy is different depending on the position of a pixel 
as seen in Fig. 9, the correction angle width may differ 
according to its degree of redundancy in order to use the data 
effectively. For example, in the high redundancy position the 
correction angle width is widened, and in the low position it 
is narrowed. 

[0048] As for the back projection process (step 105) , other than 
the conventional method, the reconstruction method proposed by 
the present applicant ( JP-A-2004-188163 ) can also be applied 
as the 3-dimension back projection process. This method is to 
simplify the conventional arcsin calculation by calculating the 
approximated line of the curve indicating the radiation source 
position to the channel direction position of the parallel beam 
projection data, and by using this method the time required for 
the reconstruction process will be drastically reduced compared 
to the conventional method. 

[0049] Also while the projection data using the circular orbit 
scanning .was described in the above-mentioned embodiment, it 
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also applicable to helical orbit scanning. However in the case 
of helical orbit scanning, in the same manner as when an object 
moves, streaks of artifacts are generated at a corresponding 
position because of the discontinuity of the data in the 
scanning data end point when only the filter-correcting 
' 2-dimensional back projection is implemented. So the data 
interpolation is carried out on the data obtained by the helical 
orbit to correct it to the circular orbit data as shown in Fig. 
2 (b) , and then the filtering correction process is implemented. 
Such data interpolation is carried out prior to the 
rearrangement to the parallel beam projection data (step 102) 
in the reconstruction procedure indicated in Fig. 4. That is, 
data such as a circular orbit is created first by implementing 
the interpolation of the data with regard to the feeding 
direction of the table corresponding to the projection data 
obtained by the helical orbit scanning. Then the CT image, in 
the same manner as the circular orbit data, can be obtained by 
implementing the above-mentioned rearrangement process (step 

102) , the correction process using the weighting function (step 

103) , the reconstruction filtering process (step 104) and the 
back projection process (step 105) to the interpolated 
projection data. 

[0050] Also the degree of the artifact generated in the imaging 
end by the helical orbit scanning is determined by the degree 
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of the discontinuity in an orbit of the X-ray source . This means 
that the degree of artifact varies according to the moving speed 
of the object (table feeding speed). Therefore in order to 

determine correction angle width index e and index F of 
projection data angle using for back-projection, it is 
necessary to consider the feeding speed of the table in addition 
to the factors to consider in the case of the normal scanning 
such as the exposure amount of the X-ray , time resolution or 
noise. For example, when the movement velocity is fast, F 
should be smaller. By making F smaller the settable maximum 

value of f also becomes smaller, but £ should be set as large 
as possible in order to obtain the corrective effect. 
[0051] While above-mentioned is the description on an 
embodiment of the reconstruction method of the present 
invention using the weighting function to the parallel beam 
projection data, the present invention is also applicable to 
the reconstruction method from the fan beam projection data. 
Hereinafter, the reconstruction method from the fan beam 
projection data will be described as the second embodiment of 
the present invention. An example of the procedure of the 
reconstruction method is indicated in Fig. 4 (b) . 
[0052] In this embodiment, the bio-movement correction of the 
projection data and the correction for eliminating the image 
distortion due to redundancy are carried out without the 
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rearrangement of the obtained fan beam projection data (step 
112) . The data correction process is for obtaining weighting 
function w( 0 , y ) to apply to the fan beam based on correction 
angle width index £ and index F of projection data angle using 
for back-projection inputted from the input unit, as in the 
first embodiment. In concrete terms, it is calculated using 
N which satisfies 2 (N_1) ^ F- £ /2<2 N by the formulas below: 

[Formula 11A] 
w ( 9 ) =0 

if [ 0 < Po 7i ] 

w ( 0 ) = (P 7 7c + 0 ) wi/ ( e 7i ) 

if [P0 7C S 9 <P X 7C , £ >0] 

w ( 9 ) =0 

if [Poti ^ 9 <P x 7t , £ =0] 

w ( 9 ) =W1*V2*2/ £ 

if [Pi 7i ^ 9 <p 2 7i , e >0 , V1=0] 
w ( 9 ) = ( ( 9 -P x 7U ) * (Wl*4/ £ ) /2 7C ) +W1*V2*2/ £ 

if [Pi 7C ^ 0 <P 2 7C , £>0,V1=£0] 

w ( 0 ) =W1 

if [PiTC ^ 0 <P 2 7C , £ =0] 

w ( 0 ) =W1 

if [PiTC ^ 0 <P 2 7C , £ =0] 

• « 

w ( B ) = ( ( 6 -P 3 71 ) *W1/ ( £ Jt ) ) +W2 
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if [P 2 tc ^ 0 <P 3 n , e >0] 



w ( 6 ) =W2 



if [P 3 tc ^ 0 <P 4 7i] 



W ( 6 ) = ( (P 4 7t - 0 ) *W1/ ( £ 7C ) ) +W2 



if [P 4 tc ^ 0 <P 5 tt: , £ >0] 



w 



( 0 ) =W1*V2*2/ £ 



if [P 5 * ^ 0 <P 6 tt r £ >0 , V1=0] 



w(0) = ((P6tc-0)*(W1*4/£ )/2tc) +W1*V2*2/ £ 



w ( 0 ) =W1 



if [P 5 7C ^ 0 <P 6 7T ., £ >0 , VI ^=0] 



if [P 5 7C ^ 0 <P 6 7t , £ =0] 



w 



( 6 ) = (P 7 7C - 0 ) Wl/ ( £ 7C ) 



w( 0 )«0 



w( 0 )=0 



if [P 6 7C ^ 0 <P 7 7T , £ >0] 



if [P 6 7C ^ 0 <P 7 7T , £ =0] 



if [P 7 7i ^ 0 ] . 



Here, the respective parameters in the above formulas are 
determined by the following respective formulas : 



[Formula 11B] 



Vl= £ -F+2 



(N-l) 



if [ £ -F+2 (N " 1} >0] 



V2- £ /2-V1 
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M=2 N 
W1=1/2 N 



W2-1/2 <N_1) 



if [ e ^0] 



W2=2*W1 



if [ £ >0 , F<M] 



W2=(2* (M-F) + £ ) *W1/ £ +W1 



if [ £ >0 , M^F] 



AA=-F 



BB=-F+ £ 



CC=M-F 



DD=M-F+ £ 



EE=F-M- £ 



FF=F-M 



GG=F- £ 



HH=F 



Po=AA 



Pi=BB 



Pi=EE 



Pi=BB 



P 2 =BB 



P 2 =EE 



P 2 =CC 



P 3 =FF 



P 3 =CC 



P 3 =EE 



P 4 =CC 



if 



if 



if 



if 



if 



if 



if 



if 



if 



if 



F<M/2+ £ /2 ] 



M/2+ £ /2^F<M/2+ £ ] 



M/2+ £ ^F] 

M/2 + £ /2^F<M/2+ £ ] 

M/2 + £ ^F<M+ £ /2] 

M+ £ /2^F] 

M/2 + £ /2^F<M] 

M^F<M+ £ /2] 

M+ £ /2^F] 

M/2 + £ /2^F<M] 
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P 4 =FF 


if 


[M^F<M+ e /2] 


P 4 =DD 


if 


[M+ £ /2^F] 


P 5 =GG 


if 


[M/2+ f /2^F<M/2+ e ] 


P 5 =DD 


if 


[M/2 + € ^F<M+ e /2] 


P 5 =FF 


if 


[M+ £ /2^F] 


Pg=GG 


if 


rF<M/2+ £ /2 1 


P 6 =DD 


if 


[M/2+ £ /2^F<M/2 + £ ] 


P 6 =GG 


if 


[M/2+ £ ^F] 


P 7 =HH. 






[0053] 


Alternatively , weighting 


function Wf ( 0 , y ) 



expressed with the formula (10) by using the same mode of 

expression as weighting function Wp ( 6 ) with regard to the 
above-mentioned parallel beam, and also with the formula (10') 
using the non-linear sub weight. 



[Formulas 12] 



Wf(0,Y) = G 



r 



Ws 
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Wf(0,y) = G< 
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[0054] An example of the weighting function for the fan beam 
back projection by the present embodiment will be illustrated 
in Fig. 10. While the diagram is shown by a sinogram, by 
replacing this weighting function with the parallel beam and 

looking at it from the perspective that 7=0, the weighting 
function shape can be changed according to the value of F and 

£ as shown in Fig. 7 (a) ~ (c) and Fig. 8 (a) -(c) . Here, when 
y = 0, in the case that F = M/2+ £ /2 or F = M/2 + £ the weighting 
function shape would be as shown in Fig. 7 (c) , and in the case 

that M/2+ £ < F < M+ £ /2 it would be shown as Fig. 7 (a) . Also 
when M/2+ £ < F = M, it would be as shown in Fig. 7 (b) . Here, 
M is M = 2 N and N is 2 (N ~ 1} ^ F- £ /2<2 N . In this way , the weighting 
function can change its weighting function shape according to 

the value of F and £ , an influence caused by the movement can 
be corrected further as the value of F and £ gets bigger and 
the noise can be reduced further as the value of £ gets smaller. 
More satisfactory resulting images can be obtained by changing 

the value of £ depending on the region. 

[0055] In the next steps, the reconstruction filtering process 
is executed by superimposing the reconstruction filtering 
function in a channel direction with regard to the fan beam 
projection data being applied to the weighting function 
obtained as above (step 113) , then CT images are obtained by 
implementing the back projection process accompanying the 
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calculation of the inverse trigonometric function and the 
distance calculation between the X-ray source and the 
reconstruction pixels (step 114) . As for the reconstruction 
filtering and the back projection process, heretofore known 
methods can be adopted. 

[0056] In this embodiment also, the discontinuity in the end 
portion of the data can be corrected with a definite corrective 
effect with any value that the projection data angle using for 
back-pro j ection have . 

[0057] This embodiment may be applied either to the 
reconstruction method for performing the calculation of the 
inverse trigonometric function and the distance calculation 
between the X-ray source and the reconstruction pixels without 
executing the rearrangement process from the fan beam to the 
parallel beam as shown in Fig. 4 (b) , or to the reconstruction 
method for implementing the rearrangement process to the 
parallel beam, reconstruction filtering process and back 
projection process after correcting the fan beam projection 
data using the weighting function. Also the weighting process 
of the fan beam projection data can be used as a cutting-out 
window for the projection data in the continuous scanning or 
helical scanning. In concrete terms, it can be applied to the 

interventional radiography for performing surgery or 

* 

hyperthermia treatment enabling the use of the contrast 
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radiography or CT apparatus as a monitor. 

[0058] As for the correction angle width and the pro j ection data 
angle using for back-projection in the process using the 
above-described weighting function, the value set in advance 
as a default may be used,, or at least one of them may be set 
arbitrary by the user. Also it is possible to obtain the 
weighting function to apply in an ex-post manner relating to 
the acquired projection data, not to obtain the weighting 
function prior to the imaging based on the correction angle 
width and/or the projection data angle using for 
back-projection either being set by the user or as a default. 
In the ex-post case, the correction angle width and the 
projection data angle using for back-projection are changed and 
applied . 

[0059] It is apparent that the objective of the present 
invention is achieved from the various embodiments described 
above. However, these are intended merely for explanation and 
exemplification, to which the present invention is not limited. 

[0060] While the tomograph using the X-ray is described in the 
present embodiment, without being limited to this, the present 
invention is applicable to tomograph using neutron radiation, 
positive electron, gamma ray or light. The scanning method is 
not limited to any of the "first", "second", "third" or "fourth 
generation" systems, and it may also be applied to multi-tube 
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type CT equipped with a plurality of X-ray sources, cathode 
scanning CT , electronic beam CT or C-arm type CT . Also, the 
form of the detector can be applied to any detectors such as 
a detector being disposed on the cylinder surface with an X-ray 
source in the center, planar detector, detector disposed on the 
spherical surface with an X-ray source in the center, or 
detector disposed on the cylinder surface with the revolving 
axis in the center. 

Embodiments 

[0061] In order to compare the conventional reconstruction 
methods using a weighting function and the reconstruction 
method of the present invention, motion artifacts and the image 
SD value were obtained and evaluated. As for the geometry, the 
one of the cone beam was adopted. 

[0062] [Simulation 1] 

The image was reconstructed using the weight being created 
by fixing the index of projection data angle using for 
back-projection as F = 1 . 1 and varying correction angle width 

index s as 0.0, 0.2, 0.4, 0.6, and the generation of motion 
artifacts was verified. The weighting function shape and the 
reconstructed images are shown in Fig. 12. In the case of 

correction angle width s = 0 as shown in Fig. 12 (a) , an 
intensive motion artifact appears since the slope of the weight 
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is not maintained. In the case of £=0.2 shown in (b) which 
is comparable to the conventional weight, some improvement can 
be recognized, yet it is still insufficient. On the other hand, 
a considerable improvement of motion artifacts was recognized 
in the reconstruction method (c) and (d) of the present 
invention . 

[0063] [Simulation 2] 

The generation of motion artifacts was verified by fixing 

correction angle width index s as 0 . 4 and varying index F of 
projection data angle using for back-projection to 0.8, 0.9, 
1.0, 1.1. The generation of motion artifacts was verified also 
in the conventional reconstruction method by using the 
weighting function by which the index of projection data angle 
using for back-projection was varied in the same manner. In 
the conventional weighting function, when the index of 
projection data angle using for back-projection is determined, 
the slope of weight also is automatically determined. 

As a result, as shown in Fig. 13, intensive artifacts were 
generated where the projection data angle using for 

back-projection was in the vicinity of 2 tc in the conventional 
reconstruction method (b) , but in the reconstruction method (a) 
of the present invention the motion artifacts were stably 
reduced regardless of the projection data angle using for 
back-pro j ection . 
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[0064] [Simulation 3] 

With regard to the reconstructed image obtained in 
simulation 2, the image SD value was measured as a noise 
evaluation. As seen in the result shown in Chart 1, in the 
reconstruction method of the present invention, the noise 
amount is somewhat increased compared to the conventional 
reconstruction method in the case of the index of projection 
data angle using for back-projection being 1.0, 1.1. However, 
since the image SD and motion artifacts are in a trade-off 
relationship, as opposed to the conventional reconstruction 
method incapable of restraining motion artifacts due to 
emphasizing the image SD, it is possible to reduce motion 
artifact considerably in the reconstruction method of the 
present invention without a significant increase of noise. 



[0065] [Chart 1] 



Proj ection Data Angle 
"F" Using for 
Back-Pro j ection 


0.8 


0.9 


1.0 


1.1 


Image 
SD 

Value 


Present 
Invention 


14.3 


13.8 


12.8 


11.8 


Conventional 
Method 


14.3 


13.8 


11.3 


11.1 



Brief Description of the Drawings 

Fig. 1 is a diagram showing the general configuration of 
the X-ray CT apparatus to which the present invention is 
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applied . 

Figs. 2 (a) and (b) are diagrams explaining the imaging 
method by which the X-ray CT apparatus of the present invention 
is adopted. 

Figs. 3 (a) and (b) are diagrams showing the detector of 
the X-ray CT apparatus in the present invention, wherein (a) 
indicates a single-array detector and (b) a multi-array 
detector . 

Figs. 4 (a) and (b) are diagrams showing a procedure of 
the image reconstruction in the present invention, wherein (a) 
indicates the first embodiment and (b) the second embodiment. 

Figs. 5 (a) and (b) are diagrams illustrating the 
rearrangement from the fan beam projection data to the parallel 
beam projection data. 

Fig. 6 is a diagram illustrating the concept of the 
weighting function being adopted by the tomogram reconstruction 
method of the present invention. 

Fig. 7 (a) - (c) are diagrams showing an example of the 
weighting function for the parallel beam back projection 
adopted by the tomogram reconstruction method of the present 
invention . 

Figs. 8 (a) ~ (c) are diagrams showing an example of the 
weighting function for the parallel beam back projection 
adopted by the tomogram reconstruction method of the present 
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invention . 

Fig. 9 is a diagram illustrating the case of applying the 
weighting functions that are different in every row in a CT 
apparatus using a multi-array detector. 

Fig. 10 is a diagram showing an example of the weighting 
function for the fan beam back projection adopted in the 
tomogram reconstruction method of the present invention. 

Figs. 11 (a) and (b) are diagrams illustrating the 
configuration of the weighting function being adopted by the 
present invention . 

Fig. 12 is a diagram showing the relationship between the 
weighting function shape and the motion artifact in the 
reconstructed image . 

Figs. 13 (a) and (b) are diagrams showing the comparison 
of the reconstructed image by the conventional reconstruction 
method and the one by the reconstruction method of the present 
invention. 

Figs. 14 (a) ~ (c) are diagrams for explaining the 
redundancy in the X-ray CT . 

Fig. 15 is a diagram for illustrating the difference of 
the redundancy by the pixels in the case of implementing the 
back projection from the data of the minimum projection data 
angle using for back-projection being 7t+2 , ym. 

• ■ « 

Fig. 16 is a diagram showing an example of a sinogram 
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indicating the minimum complete data set. 

Fig. 17 is a diagram showing the conventional weighting 
function to which the different weight is given relating to the 
two triangle regions in the revolving axis direction, while the 
sinogram is divided into less than three regions including the 
two triangle regions on a virtual sinogram. 

Fig. 18 is a diagram for illustrating the relationship 
between the image noise and the imaging time. 

Fig. 19 (a) - (c) are diagrams for illustrating the imaging 
time and the movement of an object being examined. 



55 



